Spin relaxation in quantum dots due to electron exchange with leads 
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We calculate spin relaxation rates in lateral quantum dot systems due to electron exchange be- 
tween dots and leads. Using rate equations, we develop a theoretical description of the experimen- 
tally observed electric current in the spin blockade regime of double quantum dots. Single expression 
fits the entire current profile and describes the structure of both the conduction peaks and of the 
suppressed ('valley') region. Extrinsic rates calculated here have to be taken into account for ac- 
curate extraction of intrinsic relaxation rates due to the spin-orbit and hyperfine spin scattering 
mechanisms from spin blockade measurements. 



PACS numbers: 73.23.-b, 75.30.Hx 

During the last decade, considerable progress was 
made in development of devices utilizing spin degrees of 
freedom of electron systems, [ij promising manufacturing 
of functional devices. The ultimate usage of spin devices 
for quantum information processing [j] would capitalize 
on potentially weak coupling of spin systems with their 
environments. This coupling can be of two origins: (i) 
intrinsic to material; and (ii) device specific (geometric 
or extrinsic), due to interactions with leads, substrate, 
etc. Knowledge of spin relaxation and dccohcrcncc re- 
sulting from the coupling to environment is crucial for 
evaluating the limitations of the spin-based devices. 

The intrinsic mechanisms responsible for spin relax- 
ation in semiconductor devices arc the hyperfine inter- 
action with nuclear spins and spin-orbit interaction, for 
a review see Refs. [H, H Relaxation rates were experi- 
mentally analyzed for quantum dot devices on the basis 
of GaAs heterostructures. The resulting coherence and 
spin flip times can be as short as 10 — 100ns in zero mag- 
netic field, but arc much longer than microseconds in the 
field [3. For silicon dot structures such measurements 
have yet to be done; however the two intrinsic mecha- 
nisms are small in bulk Si, and these times are expected 
to be orders of magnitude longer, > lsecond.[l,Q 

The goal of this paper is to analyze recent stud- 
ies d, UH, Q of charge transport through double quan- 
tum dot systems. Such systems are proposed for spin 
spectroscopy and control of dot states. The spin block- 
ade. Fig. [1] of electric current through double quantum 
dots is one of the most striking signatures of the elec- 
tronic spin degree of freedom in nanoelectronic devices. 
We examine in detail the electric current through a dou- 
ble quantum dot system. We point out that while recent 
studies of relaxation processes in double dots were fo- 
cused on intrinsic mechanisms of spin relaxation, the cur- 
rent contains an essential contribution arising from the 
extrinsic mechanism of spin relaxation as a consequence 
of electron tunneling between dots and leads. Below we 
propose a theoretical description of this effect. 

We show that the current remains finite, and equal to 
/ = eFg/S, even in the absence of intrinsic spin relax- 
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FIG. 1: (Color online) (a) Lateral double dot system in 
the spin blockaded regime. In the initial configuration, the 
system has one electron on the left dot, (1,0). Electric 
current through the dot flows from right to left in a se- 
quence of steps, (b) Energy level representation of trans- 
port through the dot. The sequence of transitions via singlet 
states, (1,0) —> (1, l)s ^ (2,0)s — » (1,0), quickly transports 
electron from the right to left lead, whereas current through 
intermediate triplet state, (1,0) — > (1, 1)t, is blocked since 
subsequent transition to (2,0)s(t) state is forbidden by spin 
(energy) conservation. 



ation mechanisms. Parameter Fg is the extrinsic spin re- 
laxation rate due to electron tunneling to the leads. This 
rate strongly depends on the relative energy between the 
spin state in the dot, E ~ i?|,x, ^^^^ Fermi level, Ep, 
of the connected lead 



Ts = Fh[1 - 
A = E-Ef; 



/(A)] + 272(A), 



(1) 



where /(A) is Fermi function of electrons in the lead and 
72 is written for large separation from the Fermi level. 
The first term dominates when A > — T (T is the temper- 
ature) and describes the probability of emptying a given 
spin state back into the lead, after which the loading from 
the lead can start again. This term is linear in the tun- 
nelling rate between the quantum dot and its lead. 
The second contribution 72 to the spin flip transition 
rate in Eq. ([1]), commonly referred to as co-tunnelling, 
describes the mechanism similar to the Korringa relax- 
ation mechanism (lo| of localized magnetic moments in 
bulk conductors due to their coupling to the spins of itin- 
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erant electrons. Note that although this term contains 
two additional small parameters, T/|A| and Fr/IAI, for 
a deep level, i.e. A ^ — T, it vanishes only as l/|Ap, 
compared to the exponential decay, exp(— |A|/r), of the 
first term. We use Eq. ([1]) to analyze experimental data 
of Ref. and obtain an almost perfect agreement be- 
tween the measured current and the current exclusively 
due to extrinsic relaxation processes. 

Previously, the influence of leads on processes in 
dots has been addressed in connection to different phe- 
nomena, such as the Kondo spin relaxation co- 
tunneli ng Il2l transport and spin relaxation in a single 

dot Onlfia, 



transpp] 

Til . 13] and nuclear spin relaxation in the 
Coulomb blockade regime 

The spin blockade occurs because transport through a 
double dot, Fig. [I{a), can only flow via singlet states 
and is blocked if electron entering the right dot from 
the right lead forms a triplet with the spin on the left 
dot.[li The relevant energy states corresponding to 
this regime are shown in Fig. [H^b), where (nL,rti^)s,T 
denotes a two-electron state with nL{R) electrons on the 
left (right) dot in singlet (S) or triplet (T) configuration. 
States (2,0)s,T are separated by level spacing in the left 
dot, whereas smallness of interdot exchange leaves states 
(l,l)s,T nearly degenerate. Relaxation of the spin lifts 



the spin blockade. |20l. l2ll. l22l| After time Tg triplet state 



(1, 1)t can relax by spin flip into singlet (1, l)s, and thus 
the current / ~ e/rg can be used to experimentally deter- 
mine the lifetime of an electron spin in a device, c.f. [2^. 
The importance of the spin relaxation due to electron ex- 
change with the leads is indicated by significant currents 
on the boundaries of the spin blockade region that has 
been observed in both GaAs0 and silicon [3, M lateral 
double dots. Moreover, this mechanism may dominate 
spin relaxation in silicon quantum dots, where intrinsic 
mechanisms are weak 0, 0, Hj . 



The full Hamiltonian for the double dot, Ti = + 
H\ + V, consists of interacting electrons in the dots, H^, 
free electrons in the leads, H\, and the tunnelling between 
the leads and the dots, V. Here we do not specify the 
exact form of the Hamiltonian for strongly interact- 
ing electron states in the double dot system, which in 
principle can be written in terms of the creation, dj^g-i 
and annihilation, da a, operators in the left and right 
dots (a = L, R). We only assume that the lowest eigen- 
states of i?d have the following hierarchy of energies 
i?(2,o)s < -^(1,1) < ^(2,o)t' as illustrated in Fig. (Hb). The 
Hamiltonian of free electrons in the leads Hi is written 
in terms of the creation, cj^j^.^, and annihilation, Cako-, 
operators of electrons in lead a with momentum k, spin 
a and energy ^^k: = ^£.akcli^„c^aka ■ The 

a—L,R k.cr 

coupling between states in lead a to electron states in 
the dot is represented by the tunnelling Hamiltonian V, 
with the tunnelling probabilities Wak- 



v= J2 El^^'^^.cok. 

a=L,R. k.cr 



(2) 



The relevant dot states that are involved in electron 
current are single-particle states with an electron on the 
left dot, I TO) = dl^\0) , I jO) = 4x1 0)' and two- 
particle states with one electron on each dot, that can be 
in either singlet, | S") = :^(4t4x ~ 4x4t)I ^ )' '-'^ '-'^^ 
of the triplet spin configurations, \ Tq) = --i=(4|'il,j + 

4^4^)10), |r+) = 4^4^10), |r_) = 4x4xio)- 

The system can be in one of these states and the cor- 
responding probabilities satisfy the normalization: P^g + 
Pio + Ps + Pto + Pt+ + Pt^ = 1- We omit the (2, 0) 
states, since P(2.o) ~ due to fast escape to left lead. 
The rate equations for two-particle states in the dot are 



Ps = -(71 + 272)^^5 + (^To + Ao) + 272(^To + Pt+ + Pt_) ~ TPs , (3a) 

Pto = -(7i + ^72)i'To + ^7i(^TO + no) + ^72(/^s + i'T+ +PtJ, (3b) 

Pt^ = -(71 + 72)^^+ + 7i^T0 + ^72(^^5 + Pto) , (3c) 

Pt_ = -(71 + 72)^^^ + Wo + ll2{Ps +Pto)- (3d) 



For instance, Eq. (|3ap describes the change of the sin- 
glet state population. Ps is reduced by the transitions to 
states I T ), I i ) (with rates 71/2), but is increased by 
the reverse transitions (with rates 71 /2). It also couples 
to all triplet states, | To ) and \T±), with a single rate. 



72/2. Additionally, we include the possibility of transi- 
tion from singlet 15) = (l,l)s to (2, 0)s between the 
two dots, with rate F, by introducing —TPs term. The 
remaining equations in system ([3]) for triplet states have 
similar structure, but without the interdot transitions. 
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FIG. 2: (Color online) Electron exchange with the lead in 
lowest tunnelling orders, (a) In the first order process, a single 
electron jumps on or off the dot, with electron energy in the 
lead equal to The de-tuning from the Fermi level is 

A = -E{i,i) — Ep. (b) The second order process represents 
double electron tunnelling through a virtual state. Energy of 
participating electrons in the lead are close to the lead Fermi 
energy and can be far from the dot energy level. 



The transition rates in Eqs.(l3|) are given by the lowest 
two orders in tunnelling, Fig. [21 and describe (a) electron 
hopping on (71) and off (71) the dot 



71 - r,/(A) , 71 = [1 - /(A)] 



(4) 



and (b) double exchange of electrons between the lead 
and the dot (72), via virtual states | t 0)j I i 0), result- 
ing in the creation of an electron-hole pair in the lead. 
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Here /(^) = 1/(1 + exp(^/r)) is the Fermi function in 
the right lead, and Ff, = {2ii/h)J\fF \ WR^k\ , where A/f 
is the density of states in the lead near the Fermi level 
(A,^ ^ Ep)- The expression for 72 formally diverges at 
^fe = A, but is apphcable for |A| large compared with 
temperature, when the contribution from w A is ex- 
ponentially suppressed ,24,]: 



72 ~ Tr 



27r 



|A| > T. 



(5) 



To conclude discussion of Eqs. ©-([ni), we note that they 
were obtained from general transition rate equations for 
diagonal elements of the full density matrix. Each such 
element corresponds to an eigenstate of the full lead-dot 
system, e.g. | i ) = | ) x | doti ) with energy e^. Transi- 
tion rates between these states are. 



27r 



VfmVrr, 



iO 



Since the environment relaxes much faster than the 
dot we take the trace over the electronic configurations 
{ej;,e/} in the leads, and define 7/^ = Tr^.^e^ T-^^Pg^p". 
as the rate of transition between the dot states | dot,; } 
and I dot / ) . Here Pg is the equilibrium density matrix 
for non- interacting electrons in the leads. 

We assume that the lead-dot tunnelling rates Fr l are 
larger than the interdot rate F. Indeed, the tunnelling 



between the dots is accompanied by the emission of a 
phonon (or another excitation) that carries away the en- 
ergy difference between i?(i.i)s and £'(2,0)3- Such cou- 
pling of electron states to phonon modes results in ad- 
ditional smallness of the rate F, which is determined 
by both the overlap between electron states in the two 
dots and the matrix elements of electron-phonon cou- 
pling. However, a microscopic derivation of F is beyond 
the scope of this paper. 

In the absence of magnetic fields, states that differ only 
by spin projections, are degenerate, and we introduce 
^0 = ^'to ~ and Pti = Pt+ = Pt^ ■ We use the nor- 
malization to remove redundant Pqj and obtain a system 
of equations for only three variables, PtojPti and Ps- 

In the limit of negligible tunnelling to the left dot (F 
0), this system of equations can be diagonalized: 



p^it) + r^p,,{t) ^ J,, , T] 
Pi=PT„-Pn, ri = F„ 

P2 = 3Ps - (Pt„ + 2PtJ, F2 = F„ J2 = 0; 
= + (Pto + 2PtJ, F3 = Fc, J3-271. 



1,2,3: 

Ji = 0; 
J2 



(6) 



The first two eigenmodes describe dynamics of spin in the 
double dot system, with the spin fiip rate Fg. The last 
eigenmode is for the total occupation of the right dot by 
one electron, with characteristic charge relaxation rate 
Fc. The two rates are. 



F, =71+272, Fe = FR[l + /(A)]. 



(7) 



Here we remark on spin relaxation in a single dot. The 
relevant states are empty (Pq) and singly occupied with 
spin up/down (Py^). One finds that Eqs. ([3]) for these 
states are modified, but the dynamics of (P| — P|) and 
(P| + Px) is given by Eq. ^ with ij ~ 1 and r/ = 3 
respectively, with unchanged Fg, Fc and the source terms. 

The current through the right dot is defined as the 
rate of charge escape from the singlet state into the left 
dot, I{t) ~ erPs{t) , and the stationary solution of rate 
equations, Eqs. ([3]), gives the steady current 



I = e 



./(A) 



2FF, 



1 + /(A) 4F, + F(3 + F,/Fc) 



(8) 



This equation describes the magnitude of current through 
a lateral double quantum dot in the regime of spin block- 
ade. We plot typical current profiles in Fig. [3] as a func- 
tion of de-tuning A. In experiments A is controlled by 
gate voltages and current / is mapped as a function of 
these voltages 0, Si . 

Equation ([8]) has particularly simple form in two im- 
portant limits, 

F,»F: /=ieF^^^, (peak), (9a) 
Ts < F : / = ieF„ (valley, /(A) « 1) . (9b) 
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FIG. 3: (Color online) Current through a double dot system 
in the spin blockaded regime. The dotted (dashed) lines show 
asymptotes due to first (second) order processes that domi- 
nate peak(valley) and are given by Eq. (|8]) with Fs replaced by 
7i (272). Panel (d) shows a fit to measured current (circles) 
along a line-cut of the spin blockade peak reported in §]. 

The first limit describes the left slope and the peak of 
the current, Fig. [31 As a consequence, the dimensionless 
current I /eT is nearly independent of system parameters. 
Its meaning is straightforward: the right dot is loaded 
by one electron with probability 2/(A)/(l -I- /(A)), as 
follows from Eqs.fH]) and and the singlet state is 1/4 
of this probability, as all two-particle states are equally 
populated when singlet-triplet relaxation is fast. The 
current is proportional to the escape rate F. 

The second limit, Eq. (|9bp . describes the spin blockade 
region. Depending on parameters it can exhibit distinct 
non-exponcntially decaying tail for | Aj/T ^ 1. see Fig.|31 
that should be easily observable. The factor 1/3 comes 
from the probability of the system to be in one of the 
triplet states. Indeed, according to the stationary solu- 
tion of Eqs.ilS]) in the Hmit |A| > T, the right dot is 
definitely occupied. The probabilities of finding the sys- 
tem in the singlet and triplet states are determined by 
the ratio of F and Fg, 

For Fg F the equilibrium is P5 = Pto.i = 1/4, as we 
discussed. On the other hand, if F ^ Fg the singlet state 
is almost empty (Pg ~ Fs/3F) since it takes time l/Fg to 
populate this state from one of the triplet states while it 
quickly empties into the left dot. The triplet states are 
all equally populated, each with probability « 1/3. 

Finally, we note that the second conduction peak ob- 
served in 0, [1] is explained by complementary hole trans- 
port in the cycle (2, 1) ^ (1, 1)(t^)s ^ (2, 0)s ^ (2, 1). 
In this cycle, the spin flip relaxation between (1, 1)t.s 
states occurs due to electron exchange between the left 
dot and the left lead via intermediate states (2, 1). 



In conclusion, we presented a model for spin relax- 
ation due to electron exchange between dots and leads. 
We used it to construct a theory of current through a 
spin blockaded double quantum dot, where spin flips re- 
sult in transitions between (1, l)s and (1, 1)t states, with 
rate Fg = 71 + 272. To estimate the resulting relaxation 
times, we neglect intrinsic relaxation rate 7intr (which in 
principle can be incorporated in Eq. ([8]) by replacement 
72 ^ 72 + 7inti ), and use this theory to extract relevant 
times from experiment on silicon double dotsQ. We ob- 
tain Tg ~ 0.1 ns in the peak of the current and Tg ~ 1 /.ts 
in the valley, see Fig. 3(d). These relaxation times are 
comparable and even shorter than those due to intrinsic 
spin-orbit and hyperfine mechanisms. The coupling be- 
tween electron states in the dots and the leads has to be 
taken into account for analysis of intrinsic spin relaxation 
mechanisms from the spin blockade measurements. 
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